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3-D KINEMATICS AND TOLERANCE VARIATION ANALYSIS 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application is related to co-pending U. S. Provisional Application for Patent Serial No. 
60/220,897, entitled: SIGMUND 3D - KINEMATICS MODULE ALGORITHM, and filed July 
5 26, 2000, which is hereby incorporated by reference as though set forth in full herein. 

TECHNICAL FIELD OF THE INVENTION 

The principles of the present invention generally relate to a kinematic modeling system, and 
more specifically, but not by way of limitation, to a dynamic modeling system including kinematic 
modeling having tolerances. 

1 0 BACKGROUND OF THE PRESENT INVENTION 

Different approaches for kinematics analysis have been used in different mechanical design 
applications. Approaches for modeling an open-loop system are much simpler and more 
straightforward than modeling a closed-loop system. In modeling an open-loop system, generally, 
there is no need to solve non-linear system simultaneous equations in such a system. For a closed- 
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loop system, however, solving non-linear system simultaneous equations is required, and, thus, an 
efficient and robust algorithm becomes important. One popular approach for modeling a closed-loop 
system is using a vector loop approach, which requires forming at least one closed-loop in a 
mechanical system. However, the vector loop approach cannot be applied without knowing whether 
5 an open-loop or closed-loop system exists. 

In performing the vector loop approach, kinematic governing equations are obtained by 
forming the proper number of closed-loops in a mechanical system. In general, a vector equation 
(which is equivalent to three scale equations) is obtained for each closed-loop under consideration. 
Determining the number of closed-loops and which loop is required for obtaining solutions depends 

10 on the configuration of the mechanical system, i.e., how joints/drivers are arranged in a loop. 
However, if a mechanism becomes complicated, the kinematic governing equations will typically 
become difficult to obtain. Thus, using the vector loop approach is not particularly efficient or robust. 
Furthermore, when integrating the vector loop approach with variation analysis, the solution will be, 
undoubtably, difficult to converge (especially with large variations) because, generally, there are no 

1 5 means available to control potential axis misalignments in the mechanisms due to variations. One 
typical example of a mechanism is a planar 4-bar linkage. It has one degree of freedom in a plane, 
but has zero degree of freedom in space. Consequently, if variations cause the linkages to lie on 
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different planes, the solution will almost never converge to a reasonable accuracy. For this reason, 
the current available vector loop approach, adopted in kinematics analysis with variation analysis 
packages, cannot actually solve the kinematic governing equations with variation conditions. 
Alternatively, it is a common to practice to estimate possible results based on kinematics constraints. 
5 Therefore, adopting another approach is required to solve the kinematic governing equations with 
variation conditions. 

On the other hand, in performing the Euler parameters approach, the governing equations are 
obtained by equating joint constraints directly through a coordinate transformation using a set of 
nonlinear equations expressed by general coordinates of two connected parts. Adopting Euler 

1 0 parameters in a transformation matrix simplifies the kinematic governing equations, and also takes into 
account the flexibility of a piece part. Because the kinematic governing equations are obtained by 
equating the joint constraints directly through the coordinate transformation, potential misalignments 
in the mechanisms due to variations can be easily controlled. Moreover, because there is generally 
no need to form any closed-loops, the Euler parameters approach can thus be applied in both open- 

1 5 loop and closed-loop systems without having to know the type of the system. Additionally, the 
kinematic governing equations obtained are independent of the configuration of the mechanical 
system, i.e., howjoints/drivers are arranged in the system. Therefore, even though a mechanism 
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becomes very complicated (i.e., having multiple loops and/or multiple joints in a piece part and/or 
multiple piece parts connected to a joint), the kinematic governing equations can still be easily 
obtained. Thus, it can be said that the Euler parameters approach is efficient and robust, and also 
proper for being integrated with variation analysis packages. 



system and apply random tolerance parameters at different steps along a trajectory to determine the 
results of the modeled mechanical system. By applying random tolerance parameters at different 
steps, the piece parts are changed throughout the trajectory and a user is unable to interactively view 
and/or obtain results for a single mechanical system. 

10 SUMMARY OF THE INVENTION 

To overcome the problems of conventional modeling techniques, integration ofkinematic and 
tolerance modeling is performed to allow interactive viewing and/or obtaining results of a dynamic 
mechanical system. One embodiment of the present invention includes a system and method for 
providing graphical and analytical dynamic modeling of a multi-dimensional mechanical mechanism. 
1 5 The system includes at least one processor for computing dynamic simulation of movement of the 
multi-dimensional mechanical mechanism, and rendering the multi-dimensional mechanical mechanism. 



5 



Presently, conventional modeling systems generate kinematic models of a dynamic mechanical 
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The dynamic simulation simultaneously performs kinematic and tolerance modeling. The system 
further includes a display for displaying the rendered multi-dimensional mechanical mechanism during 
the dynamic simulation. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 A more complete understanding of the method and apparatus of the present invention may 

be obtained by reference to the following Detailed Description when taken in conjunction with the 
accompanying Drawings wherein: 

FIGURE 1 is an exemplary systemblock diagram for performing the principles ofthe present 
invention; 

10 FIGURE 2 is an exemplary flow diagram describing operations of a multi-dimensional 

kinematic and dynamic modeling simulation having tolerance variations according to the principles of 

the present invention executed by the system of FIGURE 1; 

FIGURES 3 a and 3b (collectively FIGURE 3) is a detailed flow diagram for performing the 

operations according to FIGURE 2; and 
1 5 FIGURE 4 is an exemplary graphical user interface (GUI) for building, displaying, and 

simulating a multi-dimensional kinematic model according to FIGURE 3. 
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DETAILED DESCRIPTION OF THE DRAWINGS 

The principles of the present invention will now be described more fully hereinafter with 
reference to the accompanying drawings, in which exemplary embodiments of the present invention 
are shown. This invention may, however, be embodied in many different forms and should not be 
5 construed as limited to the embodiments set forth herein; rather, these embodiments are provided so 
that this disclosure will be thorough and complete, and will fully convey the scope of the invention to 
those skilled in the art. 

Multi-dimensional kinematic and dynamic modeling is an important aspect to mechanical 
engineering design and development projects. Kinematic and dynamic modeling provides a design 

1 0 engineer invaluable information for designing a robust mechanism. By modeling a mechanism and 
performing a dynamic simulation, the design engineer may determine structural weaknesses, accuracy 
of specification, and tolerance robustness, for example. A dynamic modeling system that provides 
results for the design engineer to comprehensively determine performance of the mechanism may 
include a process for simultaneously simulating kinematic and tolerance variation analysis. The 

1 5 simulation may include using a Monte Carlo random variable technique in varying either or both of 
the kinematics and tolerances. 
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FIGURE 1 is an exemplary system block diagram 1 00 for performing the principles of the 
present invention. A computing system 1 05a includes a processor 1 1 0 coupled to a memory 115. 
The processor 1 10 is further coupled to an input/output (I/O) device 120, which may further be 
coupled to a storage medium 125. The processor 1 10 may be a single processor or multiple 
processors operating in parallel. A control device 1 27, such as a computer mouse, keyboard, and/or 
touch pad, may be coupled to the computing system 105a for controlling a software program 1 29 
operating on the processor 110. Adisplay 130 is connected to the I/O device 120 in the computing 
system 105a for displaying information generated by the processor. 

In operation, according to the principles of the present invention, the computing system 1 05a 
is utilized to generate and display a kinematic model of a mechanical component 1 3 5 . To generate 
the kinematic model of the mechanical component 135, the processor 1 10 executes the software 
program 129 stored in the memory 1 1 5 as loaded from the storage medium 125. A user of the 
computing system 1 05a controls the execution of the software program 1 29 by using the control 
device 1 27. The kinematic model of the mechanical component 1 3 5 may be generated by selecting 
and grouping, via the control device 1 27, piece parts as rendered on the display 130. The model may 
define orientation and interaction of the mechanical component 135. Alternatively, the mechanical 
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component 1 3 5 may be generated by other means as known in the art, such as using a keyboard and 
alphanumeric expressions (e.g., a software program). 

The computing system 105a maybe coupled to a network 140 whereby data packets 145 
may be communicated across the network 140 to a remote computing system 105b. The network 
5 1 40 may be a local area network (LAN), wide area network (WAN), or the Internet, for example. 
The computing system 1 05a may be a server for executing the software program 129 as accessed 
and controlled by the remote computing system 105b. 

FIGURE 2 is an exemplary flow diagram describing operations of a multi-dimensional 
kinematic and dynamic modeling simulation according to the principles of the present invention 
1 0 executed by the system 1 00 ofFIGURE 1 . The process starts at step 205 . At step 2 1 0, a kinematic 
model is defined to form a mechanical mechanism (e.g. , automobile engine valve). In defining the 
kinematic model, piece part components may be arranged or connected using predefined software 
tools, such as graphical or programming tools. For example, a graphical user interface (GUI) may 
have predefined piece part tools for rods, gears, hinges, etc. , and the piece parts may be coupled in 
15 a graphical manner. 

At step 215, reference frames (e.g., multi-dimensional coordinate systems) may be applied 
to the kinematic model. The reference frames generally include an origin or base frame that operates 
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as an absolute reference and reference frames at joints to operated as relative references. Similar to 
the defining of the kinematic model at step 2 1 0, graphical or programming tools may be used to apply 
the reference frames to the kinematic model. The reference frames are generally applied at 
connection points (i.e., joints) of the piece parts that are used to form the kinematic model. By 
5 defining the reference frames, a simulator may be utilized to operate the kinematic model having a 
^ knowledge of the relative and absolute positions of the piece parts. 

\1 At step 220, tolerance parameters are applied to the kinematic model. The tolerance 

■- i 

W parameters generally take the form of +/- X, where X is a decimal number. A technique known in 

-Q 

~ r " the art as geometric dimensioning and tolerancing (GD&T) is a generic type of generating tolerance 

■SJ 1 0 parameters and may be used to form the tolerances. The tolerance parameters may be aligned with 
J^J the reference frames located at the joints of the piece parts. Alternatively, the tolerance parameters 

may be defined with different reference frames having different origins than the reference frames. The 
order of applying the tolerances (e.g., tolerance parameters before frames or frames before 
tolerances) does not affect results of the simulation. It should be understood that the tolerance 
1 5 parameters generally have the same number of multi-dimensions as the kinematic model, so that if the 
kinematic model is three-dimensional, the tolerance parameters are defined in three dimensions 
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At st ep 22 5 , the kinematic model including the tolerance parameters may be dynamically and 
simultaneously simulated. In performing the simulation, a Monte Carlo random model simulation 
technique may be employed. The Monte Carlo simulation technique generally applies a random or 
pseudo-random function to a selected variable. For example, the simulation may select one tolerance 
5 parameter, in one or more dimension, to apply a pseudo-random function having a Gaussian or other 
distribution. 

During the simulation, the simulation drives the kinematic model in a predetermined movement 
(e.g., engine valve moving up and down) while also applying the pseudo-random function to the 
tolerance parameter. The simulation may drive the kinematic model through the same movement 

1 0 while applying the pseudo-random function to different tolerance parameters on an individual basis. 
Alternatively, multiple tolerance parameters may be applied a pseudo-random distribution during the 
simulation. It should be understood that other variations of the simulation may be employed to 
produce different statistical effects. However, by performing a dynamic simulation that simultaneously 
varies both the kinematic model and the tolerance parameters, a design engineer may view and obtain 

1 5 important information of the mechanical mechanism that is not po ssible by dynamically varying the 
kinematic model and tolerance parameters separately (i.e., serially). The process ends at step 230. 
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FIGURES 3 a and 3b (collectively FIGURE 3) is a detailed flow diagram 300 for performing 
the operations according to FIGURE 2 . The process starts at step 3 02 . A number of general steps 
are shown in dashed boxes. At step 304, individual piece parts may be received as designed by a 
design engineer using a computer aided design (CAD) software program. Alternatively, the principles 
5 of the present invention may be incorporated into a CAD system. At step 3 06, reference frames, 
including joint reference frames, coordinate systems, and datum reference frames, maybe applied 
to the mechanical mechanism may be applied by selecting datum points on the piece parts. 
Additionally, tolerance parameters may be applied to the mechanical mechanism by applying the 
tolerance parameters to the piece parts. Using the reference frames, the piece parts may be 
1 0 assembled into a final mechanism at step 308. At step 3 1 0, a mechanical mechanism is created from 
the piece parts. 

As understood in the art, the number of joint constraints for each joint depends on the joint 
type. For example, a weldjoint has sixjoint constraints. Revolution, roller, cam or prismatic joints 
have five joint constraints. Cylindrical or universal joints has four joint constraints. Spherical or 
1 5 planar joints have three joint constraints. Finally, a bearing joint has two joint constraints. 

At step 312, the simulation is initialized and commenced. The initialization may include 
initializing a Jacobian matrix (Jx) based on the number of joint constraints, drivers, and piece parts 
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involved. Drivers represent a power source that determines the motion of the mechanism and are 
added to joint axes to set the mechanism in motion in a specific way. At step 314, the Monte Carlo 
simulation is executed. The simulation may include performing a Monte Carlo simulation for individual 
or multiple tolerance variables simultaneously while performing a dynamic movement of the kinematic 
5 model for the mechanical mechanism. The simulation may end upon fault condition being met, the 
final position has been reached, or the simulation duration pre-specified has ended. 

Specifically, the methodology of the present invention is shown with regard to the particular 
steps shown in solid blocks. In general step 304, the component information is read into separate 
part windows at step 3 1 6 from the CAD software. The component information may be imported in 

1 0 the form of International Graphics Exchange Syntax (IGES) files into a operating GUI. At step 3 1 8, 
materials for all piece parts may be defined to determine the properties of parts from a database that 
store the piece part data. At step 320, mass properties including center of mass may be computed. 
For each component in the assembly, the mass properties computed are mass, volume, surface area, 
and moment of inertia, and center of mass, for example. At step 322, geometrical features for 

1 5 emulation are selected from CAD or part model points. Specifically, the geometrical features on the 
piece part may be emulated by either picking vertices on CAD data, by entering the coordinates of 
the model points, or by selecting features in the CAD data. 
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In general step 3 06, reference frames or coordinate systems are then constructed from the 
selected datum points at step 324, and valid datum reference frames designated as joint reference 
frames by selecting the model points created in step 322, if applicable, at step 326. Where 
appropriate, datum reference frames created in step324 may also be used as joint reference frames. 

At step 328, deviations and tolerances on respective parts are then defined along with 
additional geometry, e.g., vectors, needed to define tolerances and measurements at step 330. To 
simulate variation of a feature, additional geometry required, such as a line, vector, or axis is also 
created at step 330. Further, to define a measurement, additional geometry required, such as a line, 
vector, or axis may also be created at step 330. In general step 308, assembly coordinate systems 
are then assigned on the parts or subassemblies at step 332. The datum reference frames or 
coordinate systems defined in step 322 may be designated as assembly coordinate systems in the 
parts or subassemblies. The datum reference frames or coordinate systems may also be used to 
designate the joint type between the parts/sub-assemblies resulting in an assembly. 

At step 334, a determination may be made whether any assembly tolerances are to be 
applied. If yes (i.e., assembly tolerances are to be applied), then deviations or tolerance on 
components are defined with respect to datum reference frames or coordinate systems at step 336; 
otherwise, where applicable, drivers are defined at step 338. At step 340, magnitude, point or area 
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of application, and direction of force applied at the assembly may be defined. At step 342, 
measurements on subassemblies and final assembly are defined to evaluate pre-determined build 
objectives. 

At step 344, a determination is made whether an initial assembly configuration exists, i.e. , 
5 whether the initial assembly configuration is possible based upon the information given. If not, then 
a further determination is made as to whether the assembly definition requires refinement at step 346. 
If not, then the simulation ends at step 348; otherwise control reverts back to step 324 for further 
datum and reference frame selection. 

In general step 3 1 0, if the assembly configuration exi st s, then a virtual assembly is generated, 
10 e.g., using a Monte Carlo simulation at step 350. The virtual assembly may not satisfy all joint 
constraints for a closed-loop system due to variations. The general coordinates of the parts in the 
assembly may then be used for an initial trial (i.e., x (0) ) for kinematics analysis. Details regarding the 
general coordinates are discussed after discussion of FIGURE 3B. At step 3 52, a Jacobian matrix 
defining relative positions of the piece parts is initialized. The size of the Jacobian matrix is initialized 
1 5 based on the number of joint constraints, drivers, and piece parts involved. Again, the number of joint 
constraints in each joint depends on its type. 
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In general step 3 1 2, the kinematic and dynamic simulation analysis starts at step 3 54 having 
a given time (driver positions) for the assembly configuration generated by the Monte Carlo 
simulation. For a new assembly configuration, time (driver positions) are reinitialized. At step 356, 
static analyses are performed to determine the assembly float in the joints. Results depend on the 
5 direction of reaction in the joints. 

In general step 3 1 4, requisite joint constraints are added to the Jacobian matrix at step 358: 
For each j oint, j oint constraints may be associated with proper elements of the Jacobian matrix J(x) 
and position vector F(x) based on the joint type and its associated parts. It should be noted that the 
column of the Jacobian matrix where the joint constraints to be applied depends on its connected 
10 parts. At step 360, driver constraints are added to the Jacobian matrix:. For each driver, driver 
constraints are placed to proper elements of the Jacobian matrix J(x) and position vector F(x) based 
on the driver type and its associated parts. It should be noted that the column of the Jacobian matrix 
where the driver constraints to be applied depends on its related parts. 



1 5 (except ground), flexibility constraints are placed to proper elements of the Jacobian matrix J(x) and 
position vector F(x). It should be noted that the column of the Jacobian matrix where the flexibility 
constraints to be applied also depends on the part itself. At step 364, step changes for the general 



At step 3 62, flexibility constraints of parts are added to the Jacobian matrix. For each part 
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coordinates are computed, for example, using Newton-Raphson and Gaussian elimination methods. 
The step changes of general coordinates Dx^ are computed by solving a set of linear equations (A6), 
presented below, using Gaussian elimination and backward substitution with total pivoting. The 
advantage of total pivoting allows for taking singularity into account, such as multiple solutions. 

At step 3 66, after updating the general coordinates using equation ( A7), a determination is 
made at step 368 as to whether the desired precision, allowable constraint tolerances, and flexibility 
of piece parts are met. For example, the step changes of general coordinates Dx ( *° are checked to 
determine whether desired precision is met. Also, at step 368, the constraint equations are checked 
to determine whether allowable tolerances and conditions of the flexibility of piece parts are satisfied. 

If any condition is not met, control reverts back to step 3 5 8 for re-evaluating the Jacobian 
matrix J(x) and position vector F(x) using the updated general coordinates (i.e., performed Newton- 
Raphson' s iterations). However, if allowable iterations are reached and the conditions are still not 
met, it means that no solutions exist for the current assembly configuration and driver positions. And, 
general coordinates of piece parts are restored to the previous convergent solution. Otherwise, 
if all conditions described above are met at step 3 68, a further determination is made at step 3 70 as 
to whether the kinematic/dynamic simulation is to be finished. 
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If it is determined at step 370 that the simulation is not to be finished, then the time step is 
advanced, and the kinematics and dynamics simulation analysis is repeated for the current assembly 
configuration at step 372 and control reverts back to step 354; otherwise, a further determination 
is made at step 3 74 as to whether the Monte Carol simulation is finished. If not, control reverts back 
5 to step 350; otherwise, the simulation ends at step 348. 

According to Newton-Raphson's Method for nonlinear systems, a given set of nonlinear 
equations may be defined: 

f 2 (x n x 29 —,xj = 0 
f 3 (x n x 3 ,—,xj = 0 

: (Al) 

fn(Xl.X2' m ~,Xj = 0 

The solution may be converged by the iteration procedure evolved from initial solution x (0) 
10 as follows: 

x (k) = x (k-i } _ J(x (k-i) yi F ( x (k-n j (A2 ) 
where k > 1 and the matrix J(\) is called Jacobian matrix, defined as: 
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df,(x) 


df,(x)l 




dx, 


dc 2 


dx n 




df 2 (x) 


df 2 (x) 


df 2 (x) 




dXj 


dx 2 


dx n 


(A3) 




%Jx) . 


dfjx) 





dx } dx 2 dx n 



and vectors x and F(x) are defined as follows: 



hi 



m 



x = 



Tito" 
/ 2 W 



(A4, A5) 



5 Therefore, a set of nonlinear equations will be reduced to solve a sequential of linear system 
equations. 
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Jfx*"- 1 * )Ax (k) = -F(xf k -'> ) (A6) 



Ax (k) = x (k) _ x (k-i) 



(A7) 



Gaussian elimination and backward substitution with maximal or total pivoting is then applied. 
5 Gaussian elimination and backward substitution may be performed by the following 

procedure: 



Given a set of linear equations 



E, :a 11 x,+a l2 x 2 +-- + a ln x n =b, 
E 2 : a 2lXl + a 22 x 2 + ••• + a 2n x n = b 2 

E n :a nl x t +a n2 x 2 + '- + a nn x n =b„ 



(Bl) 



An augmented matrix ,4 n then formed: 



10 



A = [A b] = 



a 21 °22 



a in a i,n+l 
a 2„ a 2,n+l 



(B2) 
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where A denotes the matrix formed by the coefficients and the entries in the (//+ 1 )-th column are 
the values of b that is = b t for each / = 1, 2, n. Provided that a n * 0, the operations 
corresponding to [Ej -i^ci j} /a u )E } )-^ (Ej ) are performed for each/ = 2, 3, wto eliminate 
the coefficient of x x in each of these rows. Follow a sequential procedure for / = 2, 3, . . n- 1 and 
5 perform the operation (Ej - {a fi fa u )E i )^>(Ej) for each j = i+ 1 , /+2, provided that a u = 

0. This eliminates the coefficient ofx, in £ y for all values of/ ; = 1,2, ... j-1 (2<j<ri). Theresulting 
augmented matrix has the form: 



A = 



a il a i2 



a in a i,n + l 



0 a 22 — fl^, 



(B3) 



where, except in the first row, the values of a .. are not expected to agree with those in the original 



10 matrix A. Thismatrixrepresentsalinearsystemwiththesamesolutionsetassystem(Bl). Since 
the equivalent linear system is triangular: 



E, : a n x,+ a l2 x 2 +■■■ + a ln x„ = a, >n+/ 
E 2 : a 22 x 2 +- + a 2n x n =a 2in+I 



(B4) 
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backward substitution can be performed. Solving the «-th equation for x n gives 



*„=%^ (B5) 



Solving the (w-l)-th equation for x nA and using x n yields 



**n-l,n+l ®"n-l,n X n (B6) 
a n-l,n-l 



and continuing this process, all variables can be obtained by 



_ a i,n+l a in X n a i,n-l X n-l "' a i,i+l X i+l 
X i ~ — 

^ foreach/ =n- J,n- 2,-~ ,2,1 



(B7) 



However, the procedure will not work if one of diagonal elements a u is zero during elimination 
procedure. Therefore, pivoting strategies are required. 
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Maximal (or total) pivoting maybe performed by using the following procedure. Maximal 
pivoting at the k-th step searches all the entries a^, for i = k, k+ 1 , . . . , n, and j = k, k+ 1 , . . . , n to find 
the entry with the largest magnitude. Both row and column interchanges are performed to bring this 
entry to the pivot position. 
5 In further discussion of the Euler parameters approach (as utilized in step 366), general 

coordinates are used for relative position of any piece part with respect to a global coordinate system 
of the whole assembly (or ground — non-moving part) is represented by transformation matrix as 
follows: 

2 ( e i + e l -°- 5 ) 2 ( e 2 e 3~ e i e 4) 2 ( e 2 e 4+ e i e 3 ) *x 
2(e 2 e 3 +e,e 4 ) 2^] + e] -O.s) 2{e 3 e 4 -e l e 2 ) t y 

2(e 2 e 4 -e 1 e 3 ) 2(e s e 4 +e 1 e 2 ) 2^] +e 2 4 -0.s) t z 

0 0 0 1 

10 where rotation matrix are composed of e u e 2 , e 3 , e A called Euler parameters and t x , t y , t z are 
translation components. These seven variables are denoted as general coordinates of the piece part . 
Here, an additional variable is introduced in rotation matrix. If the part is a rigid body, additional 

constraint, e) + e\ + e\ + e 4 = 1 , is applied. 
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For a given joint Joint constraint equations, which are expressed by general coordinates of 
two connected parts, may be derived by relating two joint reference frames that are attached to the 
connected parts, respectively. For example, a revolution j oint contains five equations to constrain 
three translation degrees of freedom (OjOj = 0) and two rotation degrees of freedom (e xj 1 e zi & 



1 e zi & eyj 1 e zi & e xi 1 e^) and two translation degrees of freedom (e xi 1 OjOj & 1 OjOj). 

Although types of joints may be different, constraint equations are similar. Therefore, a set of 

nonlinear joint constraint equations expressed by general coordinates of related parts may be easily 

obtained and so are their Jacobian components. 

In general, there are two types of driver constraint equations. One is for rotation driver (Z 
and the other is for translation driver (OjOj • e 2i ). Similarly, they can be easily derived by 

relating two joint reference frames that are attached to each connected part respectively. By using 

the Euler parameters approach, the following may be achieved: 

1 . Number of independent variables (general coordinates) and constraint equations 



1 Co). A prismatic joint contains five equations to constrain three rotation degrees of freedom (e^ 



depends only on the number of joint constraints, drivers, and parts involved. The 



independent variables are independent of the number ofloops existed in the assembly 
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and how joints located in a loop as in the vector loop approach. Therefore, the Euler 
parameters approach may be applied for any generic mechanisms. 

2. Flexibility (i.e., the reverse of stiflhess) of a piece part maybe taken into account by 
introducing Euler parameters. The more flexible the piece part, the larger an offset 

5 value is between the constraint, e x 2 + e 2 + e 2 + e 4 2 and 1. It is impossible to 

consider flexibility of a piece part in the vector loop approach. 

3. The ways that joint constraint equations are set up allow a user to control joint 
mismatch or axis misalignment. This feature is used when integrating kinematics 
analysis with dimension tolerance variation simulation. It is impossible to implement 

10 this feature in the vector loop approach. 

It should be understood that upon adopting the Euler parameters approach, a user may 
perform the following modifications to the process ofFIGURE 3 to take into account special cases 
to improve robustness of the system. 

1 . Use an inverse calculation of the Euler parameters from a given transformation matrix to 
15 take into account the possibility of a special singular case. 
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2 . Take into consideration the relationship between a local transformation matrix and a 
global transformation matrix for a piece part into considerations so that more generic assembled 
sequences may be applied. 

3. Consider equivalent angular solutions during iterations. 

4 . Modify the Gaussian elimination and backward substitution with total pivoting to take into 
account a case that may have multiple solutions. 

5 . Allow a user to have control over the piece part flexibility and the constraint tolerances. 
FIGURE 4 is an exemplary graphical user interface (GUI) 400 for defining, displaying, and 

simulating a multi-dimensional kinematic model 405 according to FIGURE 3 . The kinematic model 
405, which may be imported from a CAD system, is shown in window 41 0. Elements 4 1 5a-4 1 5h 
in window 420 may be utilized to define joint types for placing reference frames on the joints 
connecting the piece parts of the kinematic model 405. Window 425 provides for setting and 
adjusting parameters for controlling the simulation. As shown, the parameters may include: (i) start 
time, (ii) duration, (iii) increment, and (iv) analysis precision. Other parameters may also be included. 
In operation, the design engineer may execute the simulation and view and/or obtain data resulting 
from dynamically moving the kinematic model 405 while simultaneously performing tolerance 
variations. 
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The previous description is of a preferred embodiment for implementing the invention, and 
the scope of the invention should not necessarily be limited by this description. The scope of the 
present invention is instead defined by the following claims. 
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